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We study single jet production in electron- nucleus collisions e~ + Na — > J + X , 
using the 1-jettiness (ti) global event shape. It inclusively quantifies the pattern of 
radiation in the final state, gives enhanced sensitivity to soft radiation at wide angles 
from the nuclear beam and final-state jet, and facilitates the resummation of large 
Sudakov logarithms associated with the veto on additional jets. Through their effect 
on the observed pattern of radiation, 1-jettiness can be a useful probe of nuclear 
PDFs and power corrections from dynamical effects in the nuclear medium. This 
formalism allows for the standard jet shape analysis while simultaneously providing 
sensitivity to soft radiation at wide angles from the jet. We use a factorization 
framework for cross-sections differential in t\ and the transverse momentum (Pj T ) 
and rapidity (y) of the jet, in the region t\ <C Pj t . The restriction T\ <C Pj t allows 
only soft radiation between the nuclear beam and jet directions, thereby acting as 
a veto on additional jets. This region is also insensitive to the details of the jet 
algorithm, allowing for better theoretical control over resummation, while providing 
enhanced sensitivity to nuclear medium effects. We give numerical results at leading 
twist, with resummation at the next-to-next-to-leading logarithmic (NNLL) level of 
accuracy, for a variety of nuclear targets. Such studies would be ideal for the EIC 
and the LHeC proposals for a future electron-ion collider, where a range of nuclear 
targets are planned. 
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I. INTRODUCTION 

The discovery of the quark-gluon plasma (QGP) in heavy-ion collisions at RHIC and 
the LHC, has made possible, for the first time, laboratory studies of quark-gluon matter 
at the high densities and temperatures that existed only a few microseconds after the Big 
Bang. One of the key pieces of evidence in the discovery of the QGP was the observed [IHH] 
suppression of high transverse momentum hadrons or jets in heavy-ion collisions compared to 
that in proton-proton collisions. This suppression can be understood in terms of the energy 
loss [9HTE] experienced by fast-moving partons propagating through the QGP plasma, formed 
during the heavy-ion collision, before emerging as final-state hadrons or jets. Such nuclear 
medium effects also induce additional radiation, associated with the energy-loss mechanisms, 
that can alter the characteristics, such as the overall jet shape, of the observed radiation in 
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the final state. Such a medium modification of jet shape or jet quenching has been proposed 
in theory [T71 [18] and has been investigated at both RHIC and LHC [19-21J, where the 
nuclear medium effects are visualized by varying the jet shape parameters such as the jet 
cone size. 

Studying the medium modification of jet shape and jet production in cold nuclear matter 
would provide independent tests of the energy-loss mechanisms and shed light on this puzzle. 
In addition, it provides new and complementary observables to study phenomena related 
to cold nuclear matter including the shadowing, anti-shadowing, EMC, and fermi-motion 
effects that affect the properties of nuclear parton distribution functions (PDFs). This can 
provide complementary information to the analysis of jet quenching associated with the 
QGP, as well as independent tests of energy loss mechanisms. For example, one of the 
puzzling results [22] observed at RHIC was that heavy meson production had the same level 
of suppression as light meson production, even though one expects heavy quarks to be less 
likely to lose energy due to medium induced effects in the QGP. Similar studies with cold 
nuclear matter could shed light on this puzzle. 

The proposed electron-ion collider (EIC) [23H25] . aims to conduct detailed studies of 
electron-ion (e-A) collisions, at higher energies and luminosities than ever before, for a wide 
range of nuclear targets. Such a facility will be an ideal laboratory for nuclear studies includ- 
ing gathering detailed information on the momentum and spatial distributions of quarks and 
gluons in the nucleon, the correlations of these distributions with nucleon spin, low Bjorken- 
x physics and the associated gluon saturation physics, and in particular the effects of the 
nuclear environment on these properties as well as nuclear medium induced effects on the 
distributions of hadrons and jets. 

A powerful way to complement the nuclear studies mentioned above is through global 
event shape analyses that characterize the detailed properties of the radiation produced in 
e-A collisions. In particular, in the study of jet distributions, global event shapes which 
depend on the properties of radiation throughout the event, can provide complementary 
information to results based on analyses that focus on the properties of the radiation in and 
near the boundary of the jet. For example, energy loss in the nuclear medium can produce 
soft radiation at wide angles from the nuclear beam and final-state jet directions. A global 
event shape will capture this wide-angle soft radiation in addition to the radiation inside and 
outside the boundary of the jet. Comparing the distributions for such global event shapes 
for different nuclei in the e-A collisions can provide vital information on the relevant nuclear 
dynamics. 

In this paper, we propose the use of such a global event shape called 1-jettiness (ti) [25] 
to study single jet production in e-A collisions, 

e" +N A ^ J + X, (1) 

where electron scatters off a nucleus Na with atomic weight A, in the deep inelastic regime to 
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produce one final state jet (J). In such processes, one usually detects the final state electron 
to determine the virtuality of the exchanged gauge boson. For sufficiently large virtuality 
of the exchanged gauge boson, the machinery of QCD factorization [27] can be used to 
separate short-distance physics from non-perturbative effects which are absorbed into long 
distance parton correlation functions. Alternatively, one can consider jet production where 
the scattered electron is unobserved. In this case, it is the large transverse momentum of 
the jet that plays the role of the hard scale in the process. Such a process has been studied 
in the past in the context of spin-dependent observables [28J. 

In this work, we consider the process in Eq.Q with an additional constraint imposed by 
the 1-jettiness event shape T\. In particular, we are interested in the differential cross-section 

da A = d3(T (C ~ + — — J + X) (2) 
dy dPj T dri 

where Pj T and y are the transverse momentum and rapidity of the jet J, respectively. The 
event shape t% restricts the radiation between the final state jet and the nuclear beam 
directions. In the limit n — > 0, the final state jet becomes infinitely narrow and only soft 
radiation (of energy E ~ n) is allowed between the nuclear beam and jet directions. Any 
energetic radiation must be closely aligned with either the beam or jet directions. This is 
schematically illustrated in Fig. [T] We restrict ourselves to such configurations by imposing 
the phase space condition 

ri <C Pj T . (3) 

The detailed properties of the radiation illustrated in Fig. [T] will be affected by the nuclear 
target in the process. For example, for larger nuclei one typically expects enhanced hadronic 
activity between the jet and beam directions. The soft radiation between the beam and jet 
directions can be affected by jet quenching or energy loss as the jet emerges from the nuclear 
medium. This is because partons produced in the hard collisions could undergo multiple 
scattering inside the large nucleus and thus lead to induced gluon radiation [TU [29j ED] when 
passing through the nucleus to form the observed hadron or jet. While such effects can be 
studied by varying jet shape parameters, the information about soft radiation at wide angles 
from the jet is often lost. The main idea advocated in this paper is to study the properties 
of the observed radiation in Fig. [TJ quantified by distributions in the configuration space 
(ti, Pj T ,y), as a probe of nuclear dynamics. In particular, the 1-jettiness t\ global event 
shape is sensitive to soft radiation at wide angles from the jet and nuclear beam directions. 
Thus, 1-jettiness allows one to study jet shapes while simultaneously providing sensitivity 
to wide-angle soft radiation. 

For processes with N final state jets, the appropriate event shape is called iV-jettiness 
(rjv)[26j, corresponding to a generalization of T\ for iV-jet events. N-jettiness has been stud- 
ied previously in the context of implementing jet vetoes in hadron collider processes. New 
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FIG. 1: Schematic figure of the process e~ + Na — > J + X in the limit t\ <C Pj t . The restriction 
T\ <C Pj T allows only soft radiation between the beam and jet directions. The factorization 
framework for this process is schematically shown in Eqs.Q and (10). 



physics analyses typically classify data by the number of hard jets observed in the final state. 
Such jet binning is an effective way to enhance signals over background processes which are 
often accompanied by additional hard jets. Vetoing additional jets restricts the phase space 
for additional radiation, giving rise to large jet-veto Sudakov logarithms that can spoil the 
convergence of perturbation theory. The iV-jettiness framework, first introduced in Rcf. [26J, 
allows for vetoes on additional jets in an inclusive manner that facilitates resummation of 
the jet-veto logarithms. In this framework, the jet-veto logarithms correspond to Sudakov 
logarithms of the form ~ a™ \n 2n (tn / Q) , where Q denotes the hard scale in the process. 
Within this context, numerical results have been obtained for beam thrust (0-jettiness) dis- 
tributions for Drell-Yan processes [3TJ [32] and Higgs production [33] , threshold resummation 
in gauge boson production with two final-state jets [31], and the jet mass spectrum for Higgs 
production with one final-state jet [33] . 

In this work, we apply the technology of the iV-jettiness formalism, previously studied 
for physical processes in hadronic collisions, to electron-nucleus scattering. In this case, the 
1-jettiness event shape T\ for the process in Eq.Q is defined as 

2q A ■ p k 2qj ■ p k 
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(4) 



where the sum is over all final state particles (except the final state electron) with momenta 
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denoted by pk- The null four- vectors qA, qj denote reference vectors along the nuclear beam 
and jet directions respectively. The choice of Q a and Qj is not unique, so long as they are 
of the same order as the hard scale in the process. Different choices of Q a , Qj correspond 
to different definitions of t±, and lead to different geometric shapes for the beam and final- 
state jets ES]- One can also appropriately choose Qj as a function of the jet algorithm 
parameters to produce jets that look very close to the jets arising from commonly used jet 
algorithms. For example, in Ref. [35J, Qj was chosen to depend on the jet size parameter 
R to produce jets consistent with the anti-Zcr jet algorithm. Thus, varying the choices for 
Qj corresponds to performing a jet shape analysis. One the other hand, keeping Qj fixed 
while varying T\ corresponds to controlling the amount of radiation near the boundary or 
far away from the jet while keeping the jet algorithm parameters fixed. A jet shape analysis 
can also be performed by being differential in an extra jet shape parameter, such as the jet 
mass. In this way, 1-jettiness gives us the flexibility to study jet shapes while also providing 
sensitivity to soft radiation at wide angles from the jet. 

In this paper, we work with specific choices for Q a , Qj and leave a jet shape analysis 
in this context for future work. For specific choices of Q a and Qj, the reference vectors 
qA and qj can be determined experimentally by a minimization condition [37] such that 
the optimal choice for qA and qj minimizes the value of T\ in Eq.Q. Such an analysis 
is similar to that employed for finding the thrust axis for the thrust event shape in e + e~ 
colliders. Alternatively, qA can be chosen along the beam axis and qj can be determined by 
employing a standard jet algorithm. However, as we discuss below, in the region T\ <C Pj t , 
the computation of T\ is insensitive to the details of the jet algorithm, up to power suppressed 
terms (26]. This feature gives analytically simpler expressions compared to methods that 
depend on the detailed properties of the jet algorithm, allowing for an easier implementation 
of higher order corrections for increased precision. 

From the definition in Eq.Q, it becomes clear that energetic particles at wide angles 
from the beam and jet reference vectors qA and qj make the largest contributions to T\. On 
the other hand, energetic radiation closely aligned with either qA or qj and soft radiation 
make relatively small contributions to T\. Thus, the region of small t±, quantified by the 
condition in Eq.([3]), corresponds to a single narrow jet with only soft radiation between the 
beam and jet directions, as illustrated in Fig. [Tj 

In this paper, for the computation of T\ in Eq.Q, we make the choices 

q A = x a Pa, qj = Pj, Q a = x A AQ e , Qj = 2P Jt coshy, (5) 

where xa denotes the nucleus momentum fraction carried by the initial parton that enters 
the hard interaction. The value of xa can be determined from momentum conservation in 



terms of the electron energy, Pj T , and y and is given later on in Eq.(19). The reference 
vector qj is just chosen to be the total jet momentum 

1j = p j = ( p J T cosh V, Pj t » p J T sinh y) > ( 6 ) 
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corresponding to the massless four- vector constructed for each value of Pj T and y in Eq.O. 
The choices for Q a and Qj correspond to twice the energy of the initial parton entering the 
hard interaction and twice the energy of the final-state jet respectively. 

Experimentally, the jet momentum Pj can be obtained by applying any standard jet 
algorithm and demanding a single jet in the final state. Note that the only information 
from the jet algorithm used to compute 7~i, is the jet reference vector qj which only depends 
on the total jet momentum. In particular, other jet algorithm parameters such as the jet 
radius R do not enter the definition of T\. Different jet algorithms will in general give 
different results for the jet momentum Pj and correspondingly yield a different jet reference 
vector in Eq.(j6|. The result extracted for Pj T depends on which final state particles are 
grouped into the jet by the jet algorithm in question. However, by restricting to the region 
Ti <C Pj T , this jet algorithm dependence becomes power suppressed [26J. This can be 
understood by recalling that the limit n — > corresponds to an infinitely narrow jet with 
any additional wide-angle radiation being restricted to be soft (of energy E ~ ri). In this 
region of phase space, different jet algorithms will find the same narrow jet. In other words, 
different jet algorithms will give the same result for the total momentum of the jet, for 
events characterized by well-separated narrow jets. Any differences in the jet algorithms 
are associated with how they treat wide-angle soft radiation which has little impact in 
determining the total momentum of the final state jet, used for the jet reference vector in 
Eq.([6]). More quantitatively, differences in the jet algorithms can cause shifts in determining 
the jet reference vector qj that are suppressed by powers of ti/P Jt . 

In the theoretical calculation of the observable in Eq.(|2]), the jet momentum is defined as 

d af 2( lA-Pk 2qj-p k 

P ' = ^ Vte{ ^ Qj~ (?) 

where the sum is over all final state particles (except the scattered electron) with momenta 
denoted by pk- This definition of the jet momentum is closely tied with the definition of r% 
in Eq.Q. In the calculation of 7~i, all final state particles (pk) are associated with either 
the qA or qj directions as determined by the minimization condition in Eq.Q. The jet 
momentum is then defined as the sum of the particle momenta (p^) associated with the qj 
direction, selected by the theta function condition in Eq.(|7]). The transverse momentum 
Pj T and rapidity y of the jet, appearing in Eq.pl), just correspond to the magnitude of the 
transverse momentum component and the rapidity of four- momentum Pj of Eq.([7]). Note 
that in the region of small t±, the total jet momentum Pj as defined in Eq.Q, will differ 
from that obtained using a standard jet algorithm only by power corrections in ti/P Jt , 
associated with how the jet algorithm treats wide-angle soft radiation. 

As discussed earlier, different choices of Q a ,Qj in Eq.([7]) can be made to change the 
geometric properties of the jet. For example, as one changes Qj in Eq.Q, the particles that 
are grouped into the jet and will change. This property can be exploited to perform a jet 
shape based analysis. In particular, one can study the dependence of Pj T as a function of Qj 
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for a fixed value of t%. By choosing Qj as a function of a jet size parameter R [35], one can 
study the energy contained in the jet as a function of its cone size. This allows us to probe 
energy loss near the boundary of the jet while still retaining information on wide-angle soft 
radiation through the value of T\. 

The dynamics of the process in Eq.Q, in the restricted region T\ <C Pj t , is dominated by 
energetic collinear emissions (E ~ Pj T ) along the nuclear beam and final state jet directions 
and soft emissions (E ~ 7~i) in all directions. A convenient framework for such processes 
is given by the Soft- Collinear Effective Theory (SCET) [3BTH3] . which is a Lagrangian and 
operator based formulation of the soft-collinear limit of QCD. The SCET naturally separates 
the physics of the disparate scales T\ <C Pj t . A resummation of the Sudakov logarithms ~ 
a™ ln 2n (ri/Pj T ), associated with the restricted radiation or equivalently a veto on additional 
jets or hard radiation, naturally arises through solutions to the renormalization group (RG) 
equations in the SCET. For the process in Eq.(fi]), the SCET framework has a well-defined 
power counting in the small parameter A 

A 2 ~ ^- (8) 

In the region of T\ ~ Pj T , corresponding to allowing hard radiation or additional jets 
between the nuclear beam and jet directions, resummation effects are no longer important 
but power corrections can no longer be neglected. In addition, the jet algorithm dependence 
is no longer suppressed. The regions T\ <C Pj t and T\ ~ Pj T can be smoothly connected via 
a matching calculation. In this work, we only focus on the resummation region T\ <C Pj t , 
leaving the matching calculation for future work. 

Recently jH], a factorization framework based on the SCET, applicable in the region 
T\ <C Pj T , was introduced for the observable in Eq.(|2]). In that work, numerical results 
at the next-to-leading logarithmic (NLL) accuracy were derived for the case of a proton 
target and the impact of non-perturbative effects in the region n ~ ^qcd were studied. In 
this work, we extend the numerical results to include a wide range of nuclear targets and 
resummation at the next-to-next-to-leading logarithmic (NNLL) accuracy. In particular, we 
give numerical results for the nuclear targets: Proton, Carbon (C), Calcium (Ca), Iron (Fe), 
Gold (Au), and Uranium (Ur). The factorization formula for the observable in Eq.Q has 
the schematic form 

- H ®B® J <g) S, (9) 



dydPj T dr 



where H, B, J, and S denote the hard function, the nuclear beam function, the jet function, 
and the soft function respectively. The hard function captures the physics of the hard 
partonic interaction that initiates the final state jet. Similarly, the jet function describes 
the dynamics of collinear energetic radiation in the final state jet and the soft function 
describes the low energy radiation throughout the event. The beam function [31] B is a 
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FIG. 2: Cross-section differential in t\ and y with NNLL resummation for a proton target, at 
Pj T = 20 GeV and center of mass energy of 90 GeV. 

nuclear matrix element and encodes the physics of parton correlations in the initial nucleus, 
collinear radiation from the initial state, and the beam remnants. The various objects in 
Eq. Q have well defined field-theoretic definitions and correspond to the various parts shown 
schematically in Fig. [T] 

In the nuclear beam function B, one can separate the physics of perturbative collinear 
initial state radiation from the non-perturbative dynamics of the initial state nucleus by 
performing an operator product expansion (OPE). At leading order in the OPE, the beam 
function can be written as a convolution between a perturbatively calculable coefficient X 
and the standard nuclear PDF /a 

B~T®f A . (10) 

The OPE is an expansion in the Q 2 s (A)/t a , where Q S (A) is a dynamical nuclear scale and 
t a ~ T\Pj T denotes the virtuality squared of the initial state parton that enters the hard 
interaction after being taken off-shell by initial-state radiation. The physics of these pertur- 
bative collinear emissions from the incoming parton, after absorbing the non-perturbative 
collinear emissions into the PDF, is contained in the coefficient X. The dependence of the 
nuclear scale Q 2 S (A) on the atomic weight A of the nucleus is typically parameterized as 



Ql(A) ~ A a Ar 



(11) 
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where the parameter a determines the scaling of Q S (A) with the the atomic weight of the 
nucleus. Note that for the simplest case of a proton target (A = 1), the nuclear scale is 
just Q 2 S (A — 1) ~ Aq CD ~ I/^am where is the nucleon radius. The power corrections 
in Q 2 S {A) / {t\Pj t ) can allow one to extract information on higher twist parton correlations 
in the nucleus and the nuclear modification of gluon radiation. Note that the size of the 
power corrections will increase for heavier nuclear targets as determined by the scaling with 
the atomic weight in Eq.(ll). The size of these power corrections for a given nuclear target 
will also increase at smaller values of t% and Pj T . Thus, by analyzing the dependence of 
data on the A,Ti, and Pj T , one can extract information on the size and properties of the 
nuclear-dependent power corrections. These power corrections will manifest themselves as 
deviations from the leading twist results of Eqs.(|9]) and (10) that have increased effects 
for heavier nuclear targets and smaller values of T\ and Pj T . Also, note that while the jet 
algorithm dependence is suppressed in powers of T\/Pj T , the nuclear medium induced effects 
are suppressed by Q 2 s (A)/{t\Pj t ). Thus, for a fixed Pj T , by going to smaller values of T\ we 
can reduce the jet algorithm dependence while increasing the nuclear medium effects. 

In the region Aqcd "C T\ <C Pj t , where T\ is perturbative, the functions H, J,Z, and S 
are all perturb at ively calculable and are independent of the properties of the initial state 
nucleus. Thus, at leading twist, the only dependence on the nuclear target comes from the 
nuclear PDF and the observable in Eq.Q becomes a direct probe of the nuclear PDFs. 
In the region where T\ ~ Aqcd, the soft function S becomes non-perturbative. This can 
be understood by recalling that the soft function describes the dynamics of soft radiation 
with energy E ~ T\. In this region, a non-perturbative model must be employed for the 
soft function and its parameters can be extracted by a comparison with data. Note that 
since the soft function is independent of the nuclear target, it is a universal function. One 
can exploit this universality to extract the non-perturbative soft function from data for the 
proton target and then use it as a known quantity for processes with other nuclear targets. 

For the purposes of illustration, in Fig. [2] we show the differential cross-section in Eq. ^ 
as a function of T\ and the jet rapidity (y) for a proton target at P Jt = 20 GeV and a 
center of mass energy of 90 GeV. This result includes resummation of the jet- veto Sudakov 
logarithms at the NNLL level of accuracy. Studying such distributions in the configuration 
space of {ri,Pj T ,y} for a wide range of nuclear targets and center of mass energies, can 
provide detailed information on the structure and dynamics of nuclei. This paper is a first 
step towards such a program of exploring nuclear physics in exclusive jet production using 
a global event shape. 

In the rest of the paper, we give details of the formalism described in this section. In 
section |TlJ we describe the kinematics of the process in Eq.Q. We also describe the result 
for a naive tree-level parton model calculation and discuss how it will be modified by per- 
turbative and non-perturbative effects. In section III, we give details of the factorization 
formula shown schematically in Eqs.([9]) and (10), describe the framework used for the soft 
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function in the non-perturbative region, and discuss power corrections. In section IV 



we 



give numerical results and plots. We make concluding remarks in section [V} Various useful 
formulae and field-theoretic definitions are collected in the appendices at the end of this 
manuscript. The reader not interested in the technical details of the factorization and re- 
summation framework, can skip section III and go directly to section IV for the numerical 
results. 



II. KINEMATICS 

We carry out our analysis in the center of mass frame defined by the electron momentum 
and the average nucleon momentum in the nucleus. The electron and nucleus momentum, 
p e and Pa respectively, take the form 1 

P» = (p e ,Pe), P^ = A( P ° e ,-p e ), (12) 

where A is the atomic weight of the nucleus and the electron momentum satisfies the on- 
shell condition pi = so that the nucleus is also treated as a massless particle P\ = 0. We 
introduce the quantity Q e which is related to the electron energy as 

P°e = \Pe\ = f , (13) 

so that the hadronic Mandelstam invariant s is given by 

s = (p e + P A ) 2 = AQ 2 e . (14) 

We introduce the light-cone vector n\ and its conjugate n\ so that we can write the electron 
and nucleus momenta as 

PI = A^n% < = (1,0,0,1), 

P^=^n% n\ = (1,0, 0,-1). (15) 

The light-cone vectors satisfy n\ = n\ = and ha ■ ua = 2. The final state jet momentum 
(Pj) is given in Eq.([6|. Pj T = \Pj T \ and y denote the transverse momentum and rapidity of 
the jet respectively. We denote the light-cone four momentum vector along the jet direction 
and its conjugate as rij and fij respectively, such that n 2 j = n 2 j = 0, fij ■ rij = 2 and 
nj = -nj. 

1 Note that in the earlier work of Ref. [44], the framework was set up in the electron-nucleus center of 
mass frame. This differs from the center of mass frame of the electron and average nucleon momentum in 



the nucleus considered here. The frame defined by Eq.(12) is the one typically used in the experimental 
analysis, allowing for a more direct comparison. 
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III. FACTORIZATION 

In this section we give the factorization formula for the process in Eq.Q in the region 
T\ <C Pj T . This formula is derived using an effective field theory approach as described 
by the SCET. However, before going into the details of the factorization framework, for 
illustration and establishing the normalization of the cross-section, we first give the lowest 
order result for the observable in Eq.(|2]), using the parton model. The lowest order parton 
model result 

d 3 er (0) 1 
dydPjrdn = ^(rOEeJ^A^AO. (16) 

is obtained from tree-level partonic process convoluted with the nuclear PDF. o"o is the 
tree-level partonic cross-section differential in Pj T and y 

" dP jT dy Qlev £2 ' y > 

and the partonic Mandelstam variables s, i, and u take the form 

S = (p e + X A Pa) 2 = XaAQI, 

t = (x A P A - Pjf = -x A AQ e P jT e- y , 

U = ( Pe - Pjf = -Q e Pj T e y . 

(18) 

From the momentum conservation condition s + t + u = 0, the momentum fraction x A is 
given by 

e y Pj T 

XA = A(Q e -e-yP jT Y (19) 



Note that from Eqs.(19) and (18), the dependence on the atomic weight A completely cancels 



out in s, t, u. Thus, for all nuclear targets, a is independent of A and is equal to the partonic 



cross-section for the case of A = 1. In other words, for the kinematics given by Eq.(12), 
the ^-dependence is isolated into the factor ^f q / A (x A ,fi) in Eq.(16), the effective nuclear 
parton distribution per nucleon. 



As seen in Eq.(16), this lowest order parton-model calculation gives a cross-section that 
is proportional to 6(ti). This simply corresponds to the fact that at the lowest order the 
final state involves a jet made up of a single quark recoiling against the final state lepton. 
A calculation of the 1-jettiness in Eq.(|4]), for this configuration, trivially gives zero since the 
only final state particle that contributes is the quark which is exactly in the direction of 
the reference vector qj. However, this parton model result is not an adequate description 
since important non-perturbative effects already come into play once T\ ~ Aq C d- Recall 
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that the soft radiation, schematically shown in Fig. [TJ has energy E ~ T\ and will give rise 
to non-perturbative effects when n ~ ^qcd- These non-perturbative effects will smear the 
S(ti) distribution in a way that cannot be captured by the naive parton-model calculation. 
A rigorous analysis requires working within a factorization framework that properly treats 
the physics associated with the scales T\ <C Pj t and incorporates non-perturbative effects. 

As described earlier, the region n <C Pj T is dominated by configurations that correspond 
to a narrow jet with low-energy radiation between the nuclear beam and jet directions. 
The dynamics of this region of phase is dominated by collinear emissions along the jet and 
nuclear beam directions and soft emissions in all directions. The corresponding physics can 
be described by formulating the problem in terms of the SCET. The relevant degrees of 
freedom are the collinear modes along the nuclear beam and jet directions and the soft 
emissions with momentum scalings and virtuality given by 

beam-collinear : (tia • p,UA' P,P ±A ) ~ -Pj t (A 2 , 1> A); P 2 ~ t iPj T i 
jet-collinear : (rij -p,nj ■ p,p ±J ) ~ -P/ T (A 2 , 1, A); p 2 ~ tiP Jt , 
soft : (n A -p,n A -p,P ±A ) ~ -P/ T (A 2 , A 2 , A 2 ); p 2 ~ r 2 , 

(20) 

where p denotes a generic four-momentum vector, the power counting parameter A 2 ~ 
Ti/Pj T , and p ±A ,p ±J denote momentum components perpendicular to the beam and jet 
directions respectively. The beam-collinear modes describe the dynamics of physics along 
the beam direction, including the beam remnants. Similarly, the jet-collinear modes describe 
the dynamics of the final state jet. The typical virtuality p 2 ~ T\Pj T of the beam and jet 
collinear modes is roughly the order of the invariant mass of the final state beam and jet 
respectively. The soft modes describe the dynamics of soft radiation of virtuality p 2 ~ r\ 
that is present throughout the event. In the region n ~ ^-qcd, the soft radiation becomes 
non-perturbative. Even at small perturbative values of T\ where T\ <C Pj t , calculations 
in fixed order perturbation theory are not reliable due to the presence of large Sudakov 
logarithms ~ a" ln 2n (n/P/J that can spoil the convergence of perturbation theory. 

A resummation of large logs and an incorporation of non-perturbative effects can be 
accomplished via a factorization framework in the SCET. This was recently done, for the 
observable under consideration, in Ref. [44] . Using the by now standard techniques in the 



SCET, the factorization formula for the for the kinematics of Eq.(|12|), is given by 
d 3 a (Tq 



dydPjTdTi A 



^2 e l f dx j ds J j dt a 



e y Pj T 



xH(xAQ e Pj T e y , n; fi H )S[x - - e -vPj ) ] 

xJ q (sj,fi]fij)B q (x,t a ,fi;fi B ) (21) 
c / t a Sj 
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where the nuclear quark beam function (B q ), up to power corrections, is given in terms of 
the nuclear PDF (fi/A) as [31J 

dz ■ fx \ 
B g (x,t a ,fi;fi B ) = / —l qi {-,t a ,fi]fi B jf i/A (z,fi B ), (22) 
J x z \z / 

where the X q% are a perturbatively calculable matching coefficients and the index i runs over 
the initial parton species in the nucleus. The one-loop matching of the quark (and gluon) 
beam functions were computed in Refs. j3U [331 SHI EH] and are given in appendix [B] Note 
that the argument of the hard function is independent of A, since the A-dependence cancels 
out in the combination xA when x is evaluated at its value determined by the delta function 



in Eq.(21). The soft function appearing in Eq.(21) is defined in terms of the generalized 



hemisphere soft function (36] as 

S (n, n; fi s ) = I dk a I dkj 8{t x - k a - kj) S(k a ,kj, /i; (jt s ). 



(23) 

The generalized hemisphere soft function S(k a , kj,/i; fis), appearing on the RHS above, is a 
function of two kinematic arguments k a , kj, corresponding to the contribution to T\ of soft 
radiation grouped with the nuclear beam and jet directions respectively, as determined by 
the 1-jettiness algorithm used to calculate T\ in Eq.Q. It is also known |36j at the one-loop 
level in fixed-order perturbation theory. 

Eqs.(21) and (22) are detailed versions of the schematic formulae in Eqs.([9]) and (10) 



respectively. The intuitive role of the hard (H), nuclear beam (B q ), jet (J 9 ), and soft (S) 
functions were discussed in section |TJ All of these objects have well-defined field-theoretic 
definitions. These definitions are given in appendix [A] for completeness. Furthermore, the 
functions H, X qt , J q , and S are independent of the nuclear target and this universality can be 



exploited in nuclear studies. The argument sj of the jet function in Eq.(21) is a measure of 
the virtuality of the parton initiating the final state jet. Similarly, the argument of the beam 
function t a is a measure of the virtuality of the initial parton entering the hard scattering. 



Eq.(22), describes the process by which the initial state parton goes off-shell by an amount 



P 2 ~ t a via initial state radiation (along the nuclear beam direction) which shifts the initial 



momentum fraction from z to x as seen in Eq.(22). The perturbative coefficient X 9 * captures 



the physics of the perturbative initial state radiation. The convolution structure between 



the hard, beam, jet, and soft functions in Eq.(21) captures the dynamics of the interplay 
between the soft-collinear factored sectors. 

The hard, beam, jet, and soft scales /jlh, /Ab> A*j, and respectively are of typical size 



~ Pj T , ~ \ij ~ ^TlPj T , [i S ~ T\. (24) 

All objects in the factorization forumla are evaluated at a common scale \x. Their evolution 



from their natural scales in Eq.(24) to the scale \i are determined by their respective renor- 



malization group (RG) equations. The RG evolution between the various scales allows for a 
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resummation of logarithms associated between the scales Pj T , t±, and Kqcd- The evolution 
equations for the various objects are given by 

H(Q 2 ,fi;fi H ) = U H (Q 2 , A*, ^h)H(Q 2 , /i H ), 
X q% (^,t a ,fx;ii B ^ = J dt' a U B (t a -t' a ,^,fi B )l qi (^,t' a ,HBj 

J q (sj,ix;fxj) = J ds'jUj(sj - s'j,fx,fxj)j q (s'j,nj), 

S(k a ,kj,n;iJ, s ) = J dk' a J dk'jU s (k a - k' a ,kj - k'j, \i, \i s ) S(k a ,kj, // s ), 

(25) 

where Uh(Q 2 , fJ>, Ho),UB(t a , /i, Ho),Uj(sj, fi, /io) and Us(k a , kj, /i, fi ) are the RG evolution 
factors, from the scale /i to the scale fi, of the hard, beam, jet, and soft functions respectively. 



The nuclear PDF f^A in Eq.(22) is evaluated at the scale /i as determined by the standard 
DGLAP evolution equations. A collection of useful formulae that determine the various RG 
evolution equations is given in appendix O 



A. Factorization in position space 



The beam, jet, and soft functions that appear in Eqs.(21), (22) and (23) depend on 



variables in momentum space. One can also rewrite the factorization formula in terms of 
position space quantities. This can often simplify its implementation since the RG evolution 



equations become multiplicative instead of the convolution structure seen in Eq.(25). The 
momentum and position space functions are related via Fourier transforms as 



l qi ( — } t a ,iJ,;iJ, B ) 



J{sj,fi;fij) 
S(k a ,kj,n;fi s ) 



2n z a 



dyj 

2tt 

dykjvk, 



4vr 2 



e^J(yj,fi;fij), 

iyk ° k ° +iVk J kj S(y ka ,y kj , f i;tis) 



(26) 



where the position-space quantities appear on the RHS above and the variables 
yta-.yJiVkaiVkj are the position space analogs of t a , sj, k a , kj respectively. Note that we 
use the same notation for a given function and its Fourier transform in order to avoid to 
much clutter in notation. A given function and its Fourier transform are distinguished by 
looking at their arguments. The corresponding position space RG evolution equations are 
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multiplicative and given by 

Z q *( — ,yta>MVB) = U B (y ta ,n,fi B )l qt (— ,y ta ,fJ. B ), 

Za Z a 

J(yj,fi;Hj) = Uj(yj,ii,iij)J(yj,iij) 
S(yk a ,ykj,n;ns) = U s (y ka ,ykj,n,vs)S(yk a ,ykj,ns), 

(27) 

where U B (y ta , \l, /i ),Uj(yj, /i, /i ), and U s (y ka ,y kj , /i, /i s ) are the Fourier transforms of 
U B {t a , [A, fio),Uj(sj, fi, /io) and Us(k a , kj, /i, fi ) respectively. The factorization formula in 
terms of position-space quantities is given by 



dydPjTdr 



x 

q,i 



A{Q e -e-vPj T y 

/dy r . y T y T 

-w^e lVTTl Uj{-^- , /I, LLj)U s {y T ,yr, fi, Vs)U B (-^-, fi, /i b ) 
27T Qj Q a 



(28) 



where we have defined 



p2 

e S l-e-yP jT /Q- (29) 
B. Non-perturbative soft function 

In the region where T\ ~ A-qcd, the soft function becomes non-perturbative since now 



Lis ~ J^qcd as seen in Eq.(24). In this region, the soft function cannot be computed 
using perturbative techniques. In this soft function model can be introduced for 

phenomenological purposes and the parameters of the model can be extracted from data. 



As seen in Eq. ( A8 ) , the field-theoretic definition of the soft function is independent of the 
nuclear target. This universality can be exploited to extract the soft function from data 
collected with a proton target and used an a known quantity in for analysis with other 
nuclear targets. 

We treat non-perturbative effects with a phenomenological model for the soft function. In 
particular, we write the momentum-space generalized hemisphere soft function, that appears 



in Eq.(23), as a convolution [301 EI] of the partonic soft function (5 part .) and a model function 



("S'mod.) aS 



S{k a , kj, Li s ) = J dk' a J dk'jS part {k 

a k a , kj k'j^s)S mo d{k' a ,k'j). (30) 
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The model function satisfies the normalization condition 



dk a dkj S mo d.(k a , kj 



1. 



(31) 



The partonic soft function <S pa rt. is simply the result of the perturbative computation of 
the soft function. The model function S mo dXk' a , k'j) is typically chosen to peak around 
k' a j ~ Aqcd, so that as expected for n Aqcd, the soft function reduces entirely to <S part . 
up to power corrections in Aqcd/t~i- This can be seen by noting that the since dominant 
contribution of S mo ± comes from its peak region k' aJ ~ Aq C d and the typical scaling of the 
soft momenta in the perturbative region is k aj j ~ T\ ^> Aqcd, an OPE of the partonic soft 
function can be performed in the limit fc a> j 3> k' a j to get 



S(k a ,kj, fis) = «5p art . (k a , k j , fi s ) + O ( QC ' P ] 

7"1 



(32) 



where the normalization condition in Eq.(31 ) was used to obtain the first term above. Thus, 



as expected, in the perturbative region t\ ^> Aqcd the model soft function of Eq.(30) reduces 



to the perturbative result «S part . and the model dependence arising through S mo d. is power 
suppressed. 



The scale dependence of the soft function in Eq.(30) is contained entirely in iS part .. There 
is no scale dependence in the model function S mo d.. Since «S part . is just the perturbative 



soft function, the convolution structure Eq.(30) correctly reproduces the perturbative scale 
dependence of the soft function. 



The soft function in position space S(y T ,y T , //), that appears in Eq.(28) and is related to 



the momentum space soft function via Eq.rt26|, is correspondingly modeled using Eq.(30) as 



S(y T ,y T ,fx s ) = J dk a J dkj J dk' a J dVje-^+W 



x S part .(&a - K, kj - k'j, Hs)<Smod.(k' a , k 



(33) 



We can further simplify by writing the momentum-space partonic soft function that appears 
above in terms position-space partonic soft function as 



<Spart.(&a _ K, kj - k'j, /j, s ) 



dytjykj j 



4vr 2 



^^-W«nA k '-"j)S piat .{y ka ,y kj ,iJis) 



(34) 



Combining Eqs.(33) and (34), the convolution in Eq.(30) becomes a simple product in po- 
sition space 



<S(yr,yr,ns) = S paTt Xyr,yT, ns)S m0 dXyr,yr), 



(35) 
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where the position-space model function S mod XyT,yr) given by 



<S m od.(yr, Vt 



dk' a / dk'je-^+^S^XKik'j). 



(36) 



Further simplifications can be made by changing the variables of integration in Eq.(36). In 
particular, we introduce new integration variables u, (, defined as 



u — k a + ^j' C — k a kj> 
to rewrite the position-space model soft function iS mo d.(?/r, y T ) as 



S mo d.(yry Vt) 



du 



,u r dc > 



(37) 



(38) 



The integration over the ( variable can be perfomed to define a new single-variable function 
F m od.( u ) an d ^s position space version F mod (y) as 



F mod .(u) 



2 <->modA 2 ' 2 



mod 



(y) 



due-^ u F mod {u). (39) 



Using Eq.(39) in Eq.(38), the position-space model soft function S(y T , y T , /is) takes the form 

S{y T ,yr,^s) = <S pSLTt Xy T ,yr, ^s)F m0 d.(yr)- (40) 



Using the above relation for the soft function that appears in Eq.(28), the factorization 
formula in terms of position-space quantities, including a parameterization of soft non- 
perturbative effects, takes the form 



d 3 a 



dydPjxdTi 



00 U H (£ 2 , /i, Hh)H(£ 2 , fl H ) 



q.l 



dx 



1 dz 



S [x 



e y P 



■It 



A(Q e -e-vPj T y A 



]-t/*m(«>^b) 



^e^Ujij^-, fi, fij)U s (y T , y T , fi, fi s )U B (^-, fi, P>b) 

^ Qj Qa 



X JJr 

z' Qc 



r^B) Spart. (Vt, Vt, Us) F mod. (Vr) 



(41) 



C. Power corrections 



The factorization formula of Eqs.(21 ) and (22 ) and its equivalent form in terms of position- 
space quantities in Eq.(28), is valid at leading order in the power counting of the SCET. Sev- 
eral types of power corrections can arise and we discuss their impact on the Ti-distributions. 
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The sizes of the power corrections, in the effective theory language, are characterized by 
ratios between the scales [Mb, /ij, Us, and Q S (A). The scalings of fin, Hb, Hj, and fis 
are given in Eq.(24). Q S (A) is a dynamical scale, often referred to as the saturation scale 
[52] . associated with multiple scatterings in the nuclear medium. It depends on the atomic 
weight (A) of the nucleus and its size is typically given by Eq.(ll), where the value of a 
determines the power law dependence. If there is no color exchange between the nucleons in 
the nucleus, a ~ 1/3 [4"oTT4"T] corresponding to the path length available for the jet parton to 
have multiple scattering in the nucleus. For the simplest case of the proton, Eq.(ll) gives 
Q 2 S (A = 1) ~ A-q CD as expected. 

A systematic analysis can be performed in the SCET to derive the operator structure of 
the various power corrections. We leave such an analysis for future works, where we will 
study in detail how the multiple scattering induced gluon radiation in the final state will 
alter the radiation pattern, in particular the T\ distribution. At the moment we discuss 
nuclear-dependent power corrections, that depend on Q S (A), and how they may be probed 
through measurements of ri-distributions. As seen in Eq.( |A7[ ), the beam function is nuclear 
matrix element and is the only source of nuclear target dependence in the factorization 
formula of Eq.(21). An operator product expansion (OPE) in Q s (A)/t a can be performed 
on the beam function where the leading term is given by a perturbative function convoluted 
with the standard PDFs as shown in Eq.(22). However, higher order terms in the OPE lead 
to the more general form of the beam function 

B*(x,t a ,nv B ) = f -W-.ta./'JA**) ft/A&^ + of®^-), (42) 

J x Z \Z J \ t a / 

where the power corrections in Q 2 s (A)/t a are associated with higher twist nuclear matrix 
elements. Recall that the beam function argument t a ~ fi B ~ t~iPj t , gives the virtuality 
of the initial parton that goes off-shell via initial state radiation before entering the hard 
interaction, as explained in section III Thus, the power correction to the beam function has 
a scaling 

QM „ AaK Qc D (43) 

t a T\Pj T 

Note that this power correction has a dependence on the nuclear atomic weight through the 
factor of A a . Thus, for heavier nuclei, the effect of these of power corrections is expected 
to be larger. A detailed study of t\ distributions over a wide range of nuclear targets can 
probe these nuclear-dependent power corrections. In particular, these power corrections will 
lead to deviations from prediction of the leading twist factorization given by Eqs.(21) and 
(22) and these deviations are expected to be larger for heavier nuclei. Also, note that the 
scaling of this power correction goes like ~ 1/(t\Pj t ) compared to the typical scaling of 
1/Pj associated with power corrections to the hard function. This corresponds to the fact 
that this power correction is probing multiple scattering or nuclear modification at the beam 
scale fi B ~ a/tiPj t . 
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Power suppressed nuclear effects can also arise from multiple hard scatterings in the 
nuclear medium. These will arise as power corrections to the hard function and after the 
soft-collinear decoupling will give rise to higher twist nuclear beam functions, which will 
then be matched onto higher twist nuclear parton correlation functions. However, these 
types of nuclear-dependent power corrections have an additional suppression of ~ 1/P 



2 

Jt- 



Thus, the dominant nuclear power corrections will arise from Eq.(42) 



From Eq.(43) we see that the nuclear-medium-induced power corrections get larger for 
smaller values of T\. On the other hand, the jet algorithm dependence is suppressed by 
powers of t\/Pj t . Thus, the 1-jettiness formalism has the advantage that in the region of 
small Ti, one can study the enhanced nuclear-medium-induced power corrections without 
much sensitivity to uncertainties typically associated with implementing the details of a jet 
algorithm. 

IV. NUMERICAL RESULTS 

In this section, we present numerical results for the differential cross-section in Eq.([2]). 
We present results for a range of nuclear targets: Proton, Carbon (C), Calcium (Ca), Iron 
(Fe), Gold (Au), and Uranium (Ur). The results are at leading order in the SCET power 
counting parameter A 2 ~ ti/P Jt and include a resummation of large logarithms in T\/Pj T up 
to the next-to-next-to-leading logarithmic (NNLL) level of accuracy, using the convention in 
Table 1 of Ref. [33J for determining the order of resummation. In the region T\ ^> Aqcd, the 
numerical results are determined entirely in terms of perturbatively calculable functions and 



the nuclear PDFs. While the nuclear-size enhanced power corrections discussed in Sec. |III C 
will be left for future work, we study in detail the nuclear modification coming from the 
leading twist nuclear PDFs. For the purpose of generating numerical results, we use the 
EPS09 nuclear PDF sets from the analysis of Ref. [S3]. We also give results in the region 
where t\ ~ ^qcd where the soft function becomes non-perturbative. In this region, we use 
a phenomenological model for the non-perturbative soft function, as described in section 



III B , and show that while different model parameter choices lead to different predictions 



in the T\ ~ Aqcd region, they all converge to the perturbative result for n ^> Aqcd 



as required. Eqs.(21) and (22), corresponding to the detailed version of the schematic 



formulae given in Eqs.(|9]) and (10) respectively, give the master factorization formula for 
the leading-twist numerical results presented in this section. Power corrections will appear 
in the data as deviations from the leading twist predictions. The scaling of such deviations 



with {A, Tx, Pj T } were discussed in section III C and are expected be larger for heavier nuclei 



Thus, in addition to probing nuclear PDFs, the leading twist numerical results presented in 
this section can serve as a baseline to probe nuclear power corrections. 
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A. Nuclear PDFs and master formula 

In order to generate numerical results, the nuclear PDFs fi/A(Z,/j), appearing in the 



factorization formula in Eqs.(21) and (22), must be modeled and extracted from data. In 
the factorization formula, no assumption is made about the form of the nuclear PDF. The 
parametric form of nuclear PDFs, their connection to nuclear structure, and the extraction 
from data is still an active area of research [531458] . In this work, we use the EPS09 nuclear 
PDFs from the analysis of Ref. |53j to generate numerical results and plots. Such an 
analysis can be repeated for different parameterizations of the nuclear PDFs and it will 
be interesting to study the resulting differences. We leave such a comparative study for 
future work and limit our analysis to only working with the PDF sets in Ref. [53J. Before 
presenting the numerical results, we describe the form of these nuclear PDFs and how they 



can be incorporated into the factorization formula in Eqs.(21) and (22). 

The momentum fraction z, appearing in the nuclear PDFs fi/A{z,fi) in Eq.(22), is such 
that at z = 1 the initial parton % carries the entire momentum of the nucleus. Typically, 
models of the nuclear PDF are such that the momentum of a parton in the nucleus does 
not exceed that of the nucleon in which it is bound. In its implementation, this corresponds 
to the assumption that the nuclear PDF falls of rapidly for z > 1/A, corresponding to the 
intuitive expectation that average nucleon momentum in the nucleus is about a factor of 1/A 
smaller than the total nucleus momentum. The simplest way to incorporate this picture is 
to view the nuclear PDF as a sum of free-nucleon PDFs in the nucleus, modified by nuclear 
correction factors. After incorporating isospin symmetry, so that the u and d quarks of the 
proton PDF are the same as the d and u quarks of the neutron PDF respectively, the nuclear 
PDFs take the form [53] 

Z A — Z 

fu/A S ° 9 (x, I*) = -7 fj) fu/p(x, //) + -T— R}{X, /i) fd/p(x, fl), 



A uy A 
Z A — Z 

fd/A S ° 9 ( x , V>) = -j R t( x , V) fd/ P (x, fj.) + A R£(x, fl) f u /p{x, n), 

f!/A Sm & /i) = Jtf M) f 9/P (x, A*), (44) 

where the fi/ p (x, //) are the standard free-proton PDFs, the Rf(x,fi) denote nuclear correc- 
tion factors arising from nuclear effects on a proton bound in the nucleus, and the nuclear 
PDFs f^^ S09 (x, fi) are defined with an overall normalization factor of 1/A to give the ef- 
fective nuclear PDF per nucleon. Note that the ff^ S09 (x, fi) vanish for x > 1 since they 
are given by linear combinations of the proton PDFs fi/ p (x,[i). The argument x in Eq.(44) 



corresponds to the parton momentum fraction of the average nucleon momentum in the 



nucleus. On the other hand, the momentum fraction z in Eq.(|22|), corresponds to the par 

fEP 

HI A 



ton momentum fraction of entire nucleus. As result, the EPS09 PDFs ff/A S09 ( z i A*) 153] are 
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related to the PDFs /^(z, /i) in the factorization formula in Eqs.(21) and (22) by 



f</A(z } p) = fF, P A S09 {Az^). 



A 



(45) 



Since (A z, fi) vanishes for A z > 1, the upper limit of the range of integration for z, 



becomes 1/A Using the relation in Eq.(45), the factorization formula in Eqs.(21) and (22) 

1 dx 



can be brought to the relatively simple form 

d 3 a 



da A (Ti,Pj T ,y) = 



dydPj T dr\ 



EPS09 



(TO 



X 



dsj I dt a 



xH(£ ,fi;fi H )J q (sj,fi]fij)l qi I —,t a ,n\HB 

V x 

tg Sj 
Qa Qj 



xS ( n - ^ - ^-,fJL;ns ) fF/A S0Q {x,n B ), 



(46) 



where the subscript EPS09 on the differential cross-section indicates that the factorization 
formula has been written in terms of the EPS09 nuclear PDFs. Note that the dependence of 
the cross-section on the nuclear target is contained entirely in the nuclear PDF ff/ A S09 {x, fx), 



as seen from Eqs.(46), (47), (29), and (17). The A-dependence completely cancels out in 



the rest of the cross-section. The cross-section in Eq.(46), for electron-nucleus scattering 



has the simple interpretation of electron-proton scattering at a center of mass energy of 
s = Ql where the proton PDF has been dressed for nuclear corrections via the replacement 



fi/p ~~ >* fT/A S ° 9 - This result is a consequence of the kinematics of Eq.(|12|), the relation in 



Eq.(|45|), and the property that ff^ S09 (A z,fi) vanishes for A z > 1. 



The lower limit of integration (x*), over the argument of the nuclear PDF fE^ {x, fig) 



in Eq.(46), is given by 



e y P 



Jt 



Q e - e~yPj T 



(47) 



The corresponding range of integration [x*, 1] is then determined by the choice of the kine- 
matic variables {Q e , Pj T , y}, defined in section |H] Thus, one can access smaller values of 



Bjorken-x by increasing Q e and decreasing Pj T and y. 

Different regions in Bjorken-x are sensitive to different types of nuclear effects. For 
example, shadowing suppresses the number density of partons in the region of small Bjorken- 
x, anti-shadowing enhances the parton density at values of Bjorken-x about ~ 0.1 , the EMC 
effect suppresses the parton density at intermediate values of Bjorken-x (> 0.2), and the 
effect from Fermi motion of the nucleons enhances the parton density at Bjorken-x values 
close to one. In Fig. |3j we show numerical results for the nuclear correction factors Rf for 
the NLO PDFs in Eq.(44) for the case of a Uranium target (A = 238, Z = 92). From these 
results, generated using the publicly available code for the EPS09 PDF set p>3J, we see that 
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FIG. 3: Nuclear correction factors i?^ /^) for the NLO nuclear PDF for a Uranium target as 



defined in Eq.(44). The subscript i runs over the parton species i = {u,d,s,g}. For the u and d 



quarks, separate R- factors are given for the valence (V) and sea quarks (S). The different curves 
in each graph correspond to different values for the scale fi. By looking at the region of small 
Bjorken-x, the different curves from the bottom to the top correspond to \x = 3 GeV (Green), 
\x = 5 GeV (Blue), \x = 10 GeV (Red), and \i = 20 GeV (Purple). These plots were generated 
using publicly available code for the EPS09 PDF set [53J. 



the shape of the Rf factors clearly show the distinct regions in Bjorken-x that are dominated 
by shadowing, anti-shadowing, the EMC-effect, and Fermi motion. These different regions 
can be probed by appropriate choices for the kinematic variables Q e , Pj T , and y to determine 



the lower limit x* of the Bjorken-x integration, as determined by Eq.(47) 



In order to illustrate the effects of the nuclear correction factors Rf(x,fi) on the cross- 
sections, in Fig. [4] we first consider the ratio of the luminosity functions, that appear in the 



tree-level cross-sections (see Eqs.(16) and (45)), for a nucleus A compared to the case of a 
proton target 



x,n) 



E q e 2 q f q % s ° 9 (x^) 



(48) 



The luminosity functions appearing in the ratio above are identical to the luminosity func- 



tions for fully-inclusive deep inelastic scattering at tree-level. Using Eq.(44) for the nuclear 



PDF fg/A (x, /■*), one can study the effect of the nuclear correction factors Rf(x,fi) on 
the ratio of the tree-level cross-section for a nucleus A compared to that of a proton target, 
through luminosity ratio in Eq.(48). In Fig. |4j we show this luminosity ratio for the case of 
a Uranium target (A = 238). We see that it captures the qualitative features of shadowing, 



24 



*Ur i(x) ° 




0.005 0.010 



0.050 0.100 



0.500 1.000 



X 



FIG. 4: Luminosity ratio for Uranium to proton using NLO PDFs for /i = 3 GeV (Green), fj, = 5 
GeV (Blue), fj, = 10 GeV (Red), and /i = 20 GeV (Purple). 

anti-shadowing, EMC effects, and Fermi motion as seen in the nuclear correction i?-factors 
in Fig. [3] However, the anti-shadowing region is completely washed out due to the isospin 
effect so that Rf; ^ 1 in the anti-shadowing region. 



B. Distributions in t\,Pj t , and y 

Once higher order perturbative effects, resummation effects, and non-perturbative effects 
from soft radiation are included, the cross-section is more complicated and is given by 



Eq.(46). In this case, a simple comparison of the the tree-level luminosity ratio in Eq.(48) is 



no longer sufficient. Instead, a comparison of the predictions from the cross-section formula 



in Eq.(46) for different nuclear targets must be carried out and is the focus of the rest of 



this section. In particular, we give numerical results for a variety of nuclear targets and 
kinematic configurations in {Q e , n, Pj T , y} and discuss their implications. 



Theoretical uncertainties to the factorization formula in Eq.(46) will arise from a trunca- 



tion of the perturbative series in the calculation of the hard (if), jet (J), beam (X), and soft 
(S) functions, higher order resummation effects not included at a given level of resumma- 
tion accuracy, and non-perturbative effects in the soft function (S). In addition, theoretical 



predictions will be affected by the standard PDF uncertainties. Corrections to Eq.(46) will 



also arise from the power corrections discussed in section HI C If one is interested in probing 



these power corrections, the uncertainties mentioned for the leading twist formula of Eq.(46) 
must be sufficiently under control. 

In order to isolate nuclear effects we will compute the ratio 



Ra{ti,Pj t ,v) 



da A (n,Pj T ,y) 



(49) 



da p (T U Pj T ,y) ' 

which compares distributions in r±, Pj T , and y for a nuclear target with atomic weight A to 
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that of a proton target. In addition to the isolation of nuclear-dependent effects, the ratio 
Ra has the advantage that many of the uncertainties in the calculation of da a and da p , as 



determined by Eq.(46), cancel in the ratio. In particular, we will show that the perturbative 
uncertainties associated with resummation and the calculation of the hard, beam, jet, and 
soft functions in fixed-order perturbation theory largely cancel in the ratio, leading to much 
smaller overall uncertainty for Ra- We also show that in the region n ~ ^qcd where 
the soft function S is non-perturbative, the dependence on the phenomenological model 
implemented to describe S largely cancels in the ratio Ra- This can be understood as a 



consequence of the fact that the soft function S in Eq.(46) is universal and independent of 
the nuclear target. 

In order to estimate the perturbative uncertainty, we employ a standard scale variation 



procedure. As seen in Eq.(46), the cross-section depends on a hard function, beam function 



coefficient, jet function, and soft function which naturally live at the scales I^b, Hj, and 



fis respectively. The typical size of these scales are given in Eq.(24). All of these objects 
are evaluated at the common scale /i using their renormalization group equations to evolve 
from their natural scales. We perform a scale variation analysis similar to that carried out 
in Ref. [32]. The nuclear PDFs are evaluated at the beam scale /xg corresponding to the 



scale at which the beam function is matched on to the nuclear PDF as shown in Eq.(22) 



or more schematically in Eq.(lO). We compute the cross-sections by choosing /i = /i# and 



make four independent choices for the relative values of the scales fin, ^b, f^j, and us 



(a) jj, 


= jj, H = r 






(6)/! 


= = V 


T 2 , 


^B = 


(c) /I 


= m = v 


T 2 , 


= 


(d) \l 


= = v 







= ry/Qan, fjLj = r>/QjTi, fX S = T Ti, 

VQan, \tj = y/QjTy, Us = r"3 ln f n, 
r"3 ln T y/Q a T U hj = y/Qjn, Us = n, 

VQan, = r"2 ln T VQm, m = n, (50) 



where £ is given in Eq.(29) and r denotes the scale variation parameter. For each of these 
choices, the scale variation parameter r is varied in the range r = {1/2,2}. One can 
estimate the perturbative uncertainty by adding in quadrature the uncertainty associated 
with the variation of each of the scales fij, and fis or by analyzing the covariance 

matrix. However, for simplicity, in this work we estimate the perturbative uncertainty as 



the envelope [32], [33] of the independent scale variations in Eq.(50). These two methods are 
expected to give similar results and a more detailed discussion can be found in [59]. 

In Fig. [5j we show numerical results for the T\ distribution for a proton target. The 



factorization formula of Eq.(46) was used to calculate this distribution for the kinematic 
configuration given by Q e = 90 GeV, Pj T = 20 GeV, and y — 0, corresponding to typical 
EIC kinematics [23]. As discussed earlier, the Ti-distribution is affected by large Sudakov 
logarithms a™ ln 2n (Yi /P/ T ) in the region t\ <^ Pj T) so that the results of fixed order pertur- 
bation theory are no longer reliable and resummation is required. These Sudakov logarithms 
are associated with the veto on additional jets, enforced by the condition T\ Pj t which 
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FIG. 5: T\ distribution for a proton target with NLL' (lower red band) and NNLL (upper green 
band) resummation for Q = 90 GeV, Pj T = 20 GeV and y = 0. A more detailed description is 
given in the text. 



restricts radiation between the hard jet and the nuclear beam direction to be soft (E ~ n), 
as shown in Fig. [T] Fig. [5] shows the result for the redistribution after a resummation of the 
jet-veto logarithms. In particular, the red (lower) and green (upper) bands correspond to 
resummation at the NLL' and NNLL level of accuracy respectively. The NLL' resummation 
corresponds to NLL resummation combined with the product of the hard, beam, jet, and 
soft functions computed at NLO and using NLO PDFs. A summary of the counting of logs 
for resummation at different levels of accuracy can be found in Table 1 of Ref. [33] • The 
red (lower) and green (upper) bands in Fig. [5j are obtained from the envelope of the scale 



variations in Eq.(50). For reference, we show solid and dashed black curves corresponding 
to the scale choices (a) in Eq.(50) for r = 1, for NNLL and NLL' resummation respectively. 
The red-dotted curve corresponds to the upper envelope of the NLL' (red) band, part of 
which is hidden by the NNLL (green) band. 

Fig. [5] shows the behavior of the cross-section as one implements jet veto by restricting 
radiation at wide angles from the final-state jet and nuclear beam directions. As T\ gets 
smaller, the final-state jet becomes narrower and wide-angle radiation becomes softer. The 
enhancement of the cross-section at small T\ is a consequence of the cross-section being 
dominated by collinear emissions along the nuclear beam and final-state jet directions and 
soft emissions everywhere. In fact, the behavior of the cross-section in fixed order pertur- 
bation theory is singular in the limit T\ — > 0. It is the resummation of the jet-veto Sudakov 
logarithms a n s ln 2n (ri/Pj T ) that tames the behavior of the cross-section at small t%. 

Hard remissions between the nuclear beam and final state jet directions are allowed for 
larger T\ ~ Pj T . These emissions are perturbartively suppressed and are not accompanied 
by large Sudakov logarithms. This part of the spectrum can be described by fixed order 



27 




r,[GeVl T^GeVl 

(a)Proton (b)C and Proton 




1.0 1.5 2.0 3.0 5.0 1.0 1.5 2.0 3.0 5.0 



nfGeVl T^GeVl 

(e)Au and Proton (f)Ur and Proton 

FIG. 6: Ti-distributions with NNLL resummation for different nuclear targets for Q e = 90 GeV, 
Pj T = 20 GeV, and y = 0. In all figures, the green (upper) band corresponds to the NNLL resumed 
result for a proton target. The lower bands in different colors are the corresponding distributions 
for different nuclear targets. 



perturbation theory. A matching calculation is required to smoothly connect the resumma- 
tion region n <C Pj T with the fixed-order perturbation theory region n ~ Pj T . We leave 
such a matching calculation for future work as the focus of this paper is on the resummation 
region. For this reason, Fig. [5] is restricted to the region of small t%. We have also not have 
shown the region n < 1 GeV, since in this region the soft function S in Eq.(46), evaluated 
at the soft scale fis ~ r i ; is affected by non-perturbative effects. We give numerical results 
for this non-perturbative soft region in section IV C| 

In Fig. [6j we show the T\ distributions with resummation at the NNLL level of accuracy 
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FIG. 7: We show ratio Ra = daA/da p in the t\ distributions for various nuclear targets compared 
to the case of the proton target. For easier visual comparison, we show the results for Carbon (C) 
and Uranium (Ur) together in subfigure (f). These results include resummation at the NNLL level 
of accuracy and are calculated at Q e = 90 GeV, Pj T = 20 GeV, and y = 0, corresponding to the 
EMC region of the nuclear PDFs. 



for a variety of nuclear targets. In all plots, the green (upper) band corresponds to the 
Ti-distribution for a proton target and the lower bands in various colors correspond to 



distributions for heavier nuclear targets. In Fig. |7J we show the ratio Ra of Eq.(49) as a 
function of T\ for various nuclei at Q e = 90 GeV, Pj T = 20 GeV, and y — 0. The scale 
variation bands in Fig. [7] are obtained by computing the ratio Ra using the same scale 



choices in da a and da p and then finding the envelope of the scale variations in Eq.(50). 



This procedure corresponds to the fact that the scales fin, ^b, ^j, Us, with typical scalings 



in Eq.(24), are determined by dynamics that are independent of the nuclear target. The 
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FIG. 8: Rapidity (y) distributions with NNLL resummation for different nuclear targets for Q e = 90 
GeV, Pj T = 20 GeV, and t\ = 1.5 GeV. In all figures, the green (upper) band corresponds to the 
NNLL resumed result for a proton target. The lower bands in different colors are the corresponding 
distributions for different nuclear targets. 



nuclear dependence only arises through the structure of the nuclear PDF which is evaluated 
at the beam scale \xb- As expected, the scale variation uncertainty is dramatically reduced 
in the ratio as seen by comparing Figs. [6] and [7} 

From Figs. [6] and [7j we see that the cross-sections for heavier nuclei are generally sup- 
pressed relative to the proton. This can be understood by noting that for Q e = 90 GeV, 
Pj T = 20 GeV, and y = 0, Eq.([47|) gives ~ 0.3 for the lower limit of integration over 



Bjorken-x in Eq.(46). From Figs. [3]and[4| we see that this corresponds to probing the nu- 
clear PDFs in the EMC region. In this region the parton density in a proton bound inside a 
nucleus is suppressed compared to that of a free proton. As seen in Fig. [7j the EMC effects 
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FIG. 9: We show ratio Ra = daA/da p in the rapidity (y) distributions for various nuclear targets 
compared to the case of the proton target. For easier visual comparison, we show the results for 
Carbon (C) and Uranium (Ur) together in subfigure (f). These results include resummation at the 
NNLL level of accuracy and are calculated at Q e = 90 GeV, Pj T = 20 GeV, and t\ = 1.5 GeV, 
corresponding to the EMC region of the nuclear PDFs. 



are larger than the perturbative uncertainty quantified by the scale variation procedure. 
Thus, the ratio Ra as a function of T\ can be a sensitive probe of such EMC effects. 

In Fig. [8] we show the rapidity distributions for various nuclear targets with NNLL re- 
summation at Q e = 90 GeV, Pj T = 20 GeV, and n = 1.5 GeV. Once again, in all figures 
the green (upper) band corresponds to the rapidity distribution for a proton target and the 
lower bands in various colors correspond to heavier nuclear targets. Here also we see the 
characteristic suppression for heavier nuclei compared to the proton target. This is shown 



more quantitatively in Fig. [9j where we show the ratio Ra in Eq.(49) as a function of ra- 
pidity for various nuclei at Q e = 90 GeV, Pj T = 20 GeV, and T\ = 1.5 GeV. The scale 
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FIG. 10: Jet transverse momentum (Pj T ) distributions with NNLL resummation for different 
nuclear targets for Q e = 90 GeV, y = GeV, and n = 1.5 GeV. In all figures, the green (upper) 
band corresponds to the NNLL resumed result for a proton target. The lower bands in different 
colors are the corresponding distributions for different nuclear targets. 



variation uncertainty is given by the width of the curves and once again we see a dramatic 
reduction of the perturbative uncertainty in the ratio Ra- The size of the suppression in 
the jet rapidity distributions for heavier nuclei, provides another measure of nuclear effects. 
As seen in Fig. [9j the deviation of Ra from unity gets larger for increasing jet rapidity (y). 



This can be understood by noting that the value of x„, as determined by Eq.(47), increases 
with the jet rapidity y. For the kinematics chosen, at y = we have ~ 0.3 and for larger 
values of y we have correspondingly > 0.3. From Figs. [3] and [4] we see that for increasing 
y, we are sensitive to the nuclear PDFs deeper into the EMC region. On the other hand, for 
more negative values of the jet rapidity y, we start becoming sensitive to the anti-shadowing 
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FIG. 11: We show the ratio Ra = doA/da p in the Pj T distributions for different nuclear targets 
relative to the proton for Q e = 90 GeV, n = 1.5 GeV, and y = 0. 



region. As seen in Fig. |4j in the anti-shadowing region, the parton luminosity function of a 
bound proton is more similar to that of a free proton. Thus, as one goes to more negative 
values of the y, one is sensitive to both the anti-shadowing and the EMC regions so that 
the net effect is a smaller suppression. The jet rapidity range y £ [—1, 1], covered in Figs. [8 
and [9j corresponds to the range x* £ [0.2,0.7]. The overall effect can be summarized by 
a decreasing Ra for increasing y, as seen in Fig. [9j Note that this is in contrast to the t\ 
distributions in Fig. [7] where Ra is relatively flat as one varies T\. This can be understood 



by noting from Eq.(47), that the value of x* is independent of n, so that we are probing 



the same regions in the nuclear PDFs for different values of T\. There is however a small 



indirect dependence on T\ through the convolution structure in Eq.(46) which can affect the 
weighting of the different regions in Bjorken-x. 
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FIG. 12: We show the ratio of Pj T distributions for Uranium relative to the proton for Q e = 
{90, 120, 140,300,800} GeV, n = 1.5 GeV, and y = 0. The different choices of Q e probe different 



ranges in Bjorken-x, as seen from Eqs.(46) and (47), yielding the different sizes and shapes for 
Ra = da a/ da p . 



In Fig. 10 we show the Pj T distributions for various nuclei at Q e = 90 GeV, y — 0, 
and T\ = 1.5 GeV. Here also we see that the cross-section is suppressed for heavier nuclei 
compared to the proton due to the EMC effects that suppress the parton density in nucleons 
that are bound inside the nucleus. In Fig. 11, we show the ratio of the Pj T distributions of 
heavier nuclei to the that of the proton. The scale variation is again dramatically reduced 
in the ratio and is given by the width of the curves. We see that the relative difference in 
the cross-sections for heavier nuclei and the proton grows with increasing Pj T . This is again 
a consequence of Eq.(47) which shows that the value of grows with P Jt . For P Jt in the 
range [10 GeV, 20 GeV], x* takes on values in the range ~ [0.1, 0.3] respectively. From Fig. Ill 
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we see that for Pj T = 10 GeV we are closer to the anti-shadowing region where the parton 
luminosity for nucleons in heavier nuclei is similar to that of a free proton. For Pj T = 20 
GeV, we are well into the EMC region where there is a significant suppression in the parton 
luminosity in heavier nuclei. As a result, we see the characteristic shape of Ra as a function 
of Pj T which indicates an increased suppression for increasing P Jt . 

In the numerical results presented so far, the kinematic configurations chosen were sen- 
sitive to the anti-shadowing and EMC regions in Figs. [3] and |4} One can also probe lower 
regions in Bjorken-x, such as the shadowing region, by choosing the appropriate kinematics. 
For illustration, in Fig. 12 we show the ratio Ra as a function of Pj T for a Uranium target 
at y = and n = 1.5 GeV for the five different values Q e = 90,120,140,300,800 GeV. 
As seen from Eq.(47), by increasing Q e , one can probe lower values of x*. For example, at 
Q e = 300 GeV and Pj T = 5 GeV we have x* ~ 0.02 which is in the shadowing region as seen 
in Figs. [3]and[4j Thus, for this kinematic choice, the integration over Bjorken-x in Eq.(46) 
covers the shadowing, anti-shadowing, and EMC regions. For Q e = 800 GeV, corresponding 
to LHeC kinematics, x* ~ 0.006 for Pj T = 5 GeV and x* ~ 0.025 for Pj T = 20 GeV so that 
might start to probe small- x saturation physics (see Ref.[60] for a recent review). In this 
case, large-x physics can be isolated by going to much larger values of Pj T . Thus, the size 
and shape of the ratio Ra as a function of Pj T and Q e can be a useful way to probe nuclear 
PDFs in different regions of Bjorken-x. Similar results can be obtained for distributions in 
the jet rapidity y and T\ as a function of Q e . 

The numerical results in Figs. [6] through 12, demonstrate that distributions in Ti,Pj T , 



and y for various nuclei and different values of Q e , can be a powerful probe of nuclear PDFs, 
complementary to measurements of structure functions in inclusive deep inelastic scattering. 
Thus, a systematic program that measures distributions of various nuclei in the configuration 
space of {Q e ,Ti, Pj T ,y} can yield detailed information about nuclear structure. 

As discussed in section [Til C these distributions will also be affected by power corrections. 
The scaling of these power corrections with the kinematic variables and their dependence 
on the nuclear targets was also discussed. In particular, the dominant nuclear-dependent 
power corrections have a kinematic scaling ~ 1/ (Yi P/ T ) rather than the typical scaling 
~ 1/Q 2 (where Q is the hard scale) in fully inclusive deep inelastic scattering. Using this 
information, deviations in the data from the leading twist predictions of Eq.(46) can be used 
as a probe of power corrections. In particular, the size of these deviations as a function of 
{A, Q e , P Jt , y, t%} can provide detailed information on the behavior and size of the power 
corrections. Such a detailed study of power corrections is left as future work. 



C. Non-perturbative soft radiation effects 



In the numerical results presented so far, we have restricted to the region T\ > 1 GeV so 
that the soft function S in Eq.(46) remains perturbatively calculable. The soft function S, 




FIG. 13: In sub-figure (a), we show the t\ -distributions for the proton and Uranium targets in the 
region t\ ~ Aqqd- The solid red (dashed blue) curves correspond to the soft function model I 
(II) in Eqs.(51) and (54). The top (bottom) two curves are for the proton (Uranium) target. In 
sub-figure (b) we show the ratio Ra = daA/da p as a function of t\ for the Uranium target using 
the soft function model I (II) as denoted by the solid red (dashed blue) curves. This plot shows 
that the model dependence of the soft function (seen in sub- figure (a)) largely cancels out in the 
ration Ra since the solid red (dashed blue) curves, corresponding to models I (II) respectively, 
largely overlap. The plots are for the kinematic configuration Q e = 90 GeV, Pj T = 20 GeV, and 
y = 0. 



which describes the dynamics of soft radiation with energy E ~ n, naturally lives at the 
scale ~ T\. Thus, in the region t\ ~ Aq C d the soft function becomes non-perturbative. 



As discussed in section III B for phenomenological purposes we implement a model for the 



soft function, as shown in Eqs.(30) and (31), as a convolution between the perturbative 



soft function 5 part . and a model function function S mo d.- As explained in section III B such 
a parameterization has the property that for T\ 3> Aqcd the soft function model reduces 



to the perturbative result <S part . as desired. From Eqs.(30) through (41), the soft function 



model can be parameterized by the function F mo d. which is related to S mo d. as in Eq.(39) 
For the purposes of generating numerical results, we employ the parameterization 



W(w) = ^ [jj Exp 



(u - bf 
A 2 



(51) 



where the values of the parameters a, b, A determine the model. The normalization N(a, b, A) 
is chosen to satisfy the condition 



duF mod (u) = 1, 



(52) 



which is equivalent to the normalization condition in Eq.(31). The parameters a, 6, A are 
chosen so that F mo d. peaks in the region u ~ Aqcd, which ensures that the soft function 
reduces to the perturbtive result for r x 3> Aq C d up to power corrections in Aq C d/ t i, as 



explained in the discussion around Eq.(32). 



36 



For the soft scale /is appearing in iS part . in the soft function model of Eq.(30), we make 
the choice 



„-min \ 2 



tis = T 1 \ll+\^-) , (53) 

with rf" n = 1 GeV. This choice has the property that in the limit T\ — > the soft scale in 
5 part . remains perturbative fig ~ * 1 GeV. For t x ^> r™ m , the soft scale reduces to \is ~ t%. 



In Fig. 13 (a), we show the Ti-distributions for a proton and Uranium target in the region 
that includes t\ ~ Aqcd at Q e = 90 GeV, Pj T = 20 GeV, and y = 0. Two curves are shown 
for the proton (top two curves) and Uranium (bottom two curves) targets. The two curves 
for each target correspond to using the two sets of model parameters 

Model I: a = 2.0, b = -0.2, A = 0.2 GeV, 

Model II: a = 1.2, b = —0.1, A = 0.3 GeV, (54) 
where Model I and Model II correspond to the solid-red and blue-dashed curves respectively 



in Fig. 13 (a). We see that for small values of T\ ~ Aqcd, there is a significant difference 
in the distributions. However, in the region t% > 1 GeV, the curves for models I and II 
converge to the perturbative result as expected. 

Since the soft function is universal and independent of the nuclear target, the model 
parameters a, b, A can be extracted from measurements of the T\ distributions in the region 
Ti < 1 GeV using a proton target. Similarly, one can also measure distributions in y and 
Pj T in the region T\ ~ ^qcd i n order to extract the soft function parameters. Once the 
parameters are extracted from data on the proton target, the soft function model can be 
used as a known input for the heavier nuclei. 



In Fig. 13 (b), we show the ratio Ra for the Uranium target as a function of t±. Again, the 



solid- red and blue-dashed curves correspond to using models I and II in Eq.(54) respectively. 
We see that in the ratio Ra, the model dependence is greatly reduced as seen by large overlap 
of the two curves corresponding to the two different model soft functions. Thus, in addition 
to the reduction in the perturbative uncertainties, there is also a reduction in the uncertainty 
associated with the non-perturbative soft function, when considering the ratio Ra- 



V. CONCLUSIONS 



In this paper, we studied electron-nucleus collisions with one final state jet e~ + Na — > 
J + X, as a probe of nuclear structure and dynamics. We used a factorization framework 
to calculate the cross-section differential in 1-jettiness (rj) and the transverse momentum 
(Pj T ) and rapidity (y) of the jet. The 1-jettiness variable T\ is a global event shape that 
quantitatively characterizes the degree to which the final state looks 1-jet-like and provides 
enhanced sensitivity to radiation at wide angles from the jet and nuclear beam directions. 
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A veto on additional jets is imposed by restricting to the region T\ <C Pj t , allowing only 
soft radiation (E ~ T\) between the beam and jet directions. This phase space restriction 
induces Sudakov jet-veto logarithms ~ ct n In 2 ™ (ti/P Jt ) that can spoil the convergence of 
perturbation theory and requires resummation. Implementing the jet-veto and performing 
the resummation using the 1-jettiness global event shape, allows for better theoretical control 
compared to methods that depend on the details of a jet algorithm. This allows one to 
perform analytic calculations at higher orders in perturbation theory and resummation, 
allowing for improved accuracy. 

Distributions in t\ provide a quantitative measure of the hadronic activity or the pattern 
of radiation between the beam and jet directions. By studying this distribution for a range 
of nuclear targets and at different kinematics, one can probe various aspects of nuclear 
physics. At leading twist, the factorization formula takes the schematic form in Eqs.Q 
and (10) and directly probes the nuclear PDFs. This can be understood by noting that in 
the region Aq C d "C T\ <C Pj T , the leading-twist cross-section is given entirely in terms of 
perturbatively calculable universal functions (independent of nuclear target) and the nuclear 
PDFs. Thus, comparing distributions in {n, Pj T , y} for a range of nuclear targets and center 
of mass energies allows for a systematic study of the nuclear PDFs. 

Power corrections beyond leading twist will probe dynamical nuclear effects such as higher 
twist correlations and nuclear modification effects such as jet quenching and energy loss 
mechanisms of fast-moving partons through cold nuclear matter. We gave a discussion 
of the various sources of power corrections and the ratios of energy scales that determine 
their sizes based on the power counting. The effective field theory framework allows one 
to systematically derive operator definitions of the power corrections to the leading twist 
factorization formula. We leave such a detailed study of power corrections for future work. 
From a phenomenological point of view, the size and shape of the various power corrections 
can be investigated by looking at the size of deviations between data and the leading twist 
prediction for a range of nuclear targets and kinematics. One typically expects that the 
nuclear-medium-induced power corrections will have a much larger effect for heavier nuclei 
and one might correspondingly expect larger deviations from the leading twist predictions 
for heavier nuclei. 

As a first step, in this paper we give numerical results at leading twist with resummation 
at the next-to-next-leading (NNLL) logarithmic order for distributions in {n, Pj T ,y} for 
the nuclear targets: Carbon, Calcium, Iron, Gold, and Uranium. We also give results for 
the ratio of these distributions between heavy nuclei and the proton. We find that there 
is a dramatic reduction in the scale variation uncertainty, as expected, when considering 
such ratios of distributions. In the region t\ ~ ^qcd, the soft function in the leading twist 
formula becomes non-perturbative since the energy of the soft radiation has the scaling 
E ~ T\. In this case, we employed a model for the soft function such that it has the correct 
renormalization group properties and reduces to the perturbative result for T\ Aqcd- 
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Furthermore, this soft function is independent of the nuclear target. This universality can 
be exploited to extract the non-perturbative soft function from data on the proton target and 
then used for the case of heavier nuclei. We also showed that in the ratio of distributions 
between a heavy nucleus and the proton, the dependence on the parameters of the non- 
perturbative soft function model largely cancel. 

Our leading twist numerical results indicate that distributions in {ti, Pj T ,y} are quite 
sensitive to differences in the nuclear PDFs. By choosing appropriate kinematics one can 
probe various regions in Bjorken-rr of the nuclear PDFs. This allows one to conduct studies 
of nuclear phenomena such as shadowing, anti-shadowing, and the EMC effect. 

We conclude by noting that this is just a first step in using event shapes for exclusive jet 
production as a probe of nuclear dynamics. There are many further directions to pursue, 
including constructing new observables that are variants of the one we studied in order to 
probe nuclear dynamics in different ways. For example, one can be separately differential 
in the contributions to 1-jettiness from the beam and jet regions, employ a standard jet 
shape analysis while still retaining information about wide angle soft radiation, construct 
analogous observables in different reference frames, study exclusive multi-jet production, 
and perform detailed studies of the various nuclear-medium-induced power corrections. One 
might also consider extending such methods to studies of p-A and A-A collisions. We look 
forward to such further developments which can be part of the broad program of physics 
envisioned by the EIC and LHeC proposals for a future electron-ion collider. 
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Appendix A: Field-theoretic definitions 

The electromagnetic quark current at the jet production vertex for flavor q is given by 

W = ^7^,(0), (Al) 
and is matched onto the an operator in the SCET as 

^(0) = du A rfwj C(wiUj,/i) v/J^M ^ A X 9 ,^(0), (A2) 
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where C(uaWj, [x) is the Wilson coefficient that contains the physics of the hard scale and 
is related to the hard function H that appears in Ea.(21) as 

H(u) A u)j,ii) = \C(u a ojj,h)\ 2 - (A3) 

The variables ua,j are label momenta that denote the large light-cone momentum compo- 
nents of the collinear quark fields along the nuclear beam and jet directions respectively. 
The fields \ denote collinear quark fields £ dressed by collinear Wilson lines W that sum up 
collinear emissions 



Xn(x) = Wlxn(x), 



fi-AZJx) 



perms. 



n ■ V 



(A4) 



and the Y denote soft Wilson lines that sum up eikonal soft emissions 



Y nA {x) = PExp 



ds ua ■ A s (s n h 



Y nj {x) =PExp 



W 



ds rij ■ A s (s n^j 



(A5) 



For detailed explanations of the notations used above, we refer the reader to the original 
SCET papers in Refs. 



The quark jet function J q in Eq.(21) is defined as 

1 



J q (sj = Wjr + ,/i) 



Im 



i / /k''- 1 <0|T{x n , Wj (0)^x nj (x)}|0> 



(A6) 



and the beam function is defined as 



B«(x,t,ii) 



1 



2xi%a ■ Pa J 47r 



db- 



e l 2,t A . VA J2(pA\XnJ{xn A -Pa -V j )(b )^Xn A (0)\PA)- 



pols. 



(A7) 



Finally, the soft function in Eq.(21) is defined in terms of the generalized hemisphere soft 
function S(k a ,kj, fi) through Eq.(23). The definition of S{k a , kj, /i) is given by 



S(k a ,kj,n) = ±^Tr(0|T[^y nj ](0)5(A: a -^^)^ 



L x fl 



Nr. 



x s 



Q. 



x (X s \T[YlY nA ](0)\0) 

where we have defined 

<J) ST^ _ nt 2( lA-Pk 2qj-Pk* 



(A8) 



K 



x s 



kex s 



Qc 



Q.. 



K 



(«) 
\ 



2qj ■ p k 2q A ■ p k . 



kex s 



Q.. 



Qa 



(A9) 
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Appendix B: Fixed order NLO expressions 



by 



In momentum space and at tree-level, the hard, jet, beam, and soft functions are given 

J^ \sj,lij) = S(sj), 
X^{^,t a ^j) = 5 qi 5(l-^)5(t a ), 



<sSrt.(k a ,kj,fi S ) = 6(k a )6(kj). 



The NLO expressions for the hard 



and jet [63, 64J functions are given by 



J q{1) (s,Li) 



O^sCf 

4:71 



-21n 2 ^- + 61n^ 



S(s)[7- ?T 



16 



71 



3 J 



\i A V s 



+ 



Li 2 I 



/x 2 fl(s) ln(s//i 2 



and the beam function coefficients [48J and the soft function [36] are given by 



2tt 



R 6{l-x)- 
b 1 — x 



In x + (1 — x) 



+ S(t) 



ln(l — x) 
1 — x 



(1 + x 2 ) + 



5{l-x) 



+ 



1 r/i 2 i l + x 2 



t 



2tt I /i 2 



(l-x) + ^ 

(1 - 2x + 2x 2 ) 



1 — x 



+ 5(t) (1 - 2ar + 2x 2 )(ln 1) + 1 



x 



(Bl) 



(B2) 



(B3) 



*5part.(^a? kj, fl) 



OL s Cp J 8 |"0(/c a )//ln fe a //i] 8 5(/c a ) ^9(kj)jHnkj/p, 



7T 



47T /i 



5(fca)5(fci) 



k n 



(B4) 



where we have defined 



fx = fJ,y/SaJ, 



SaJ 



Zqa ■ qj 

QaQj 



(B5) 
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In generating numerical results we worked with the position-space version of the factor- 
ization formula as given in Eq.(41). The jet, beam, and soft functions in momentum space 

analogs as 

e*>"J{yj,li;Hj), 



are given in terms of their position-space analogs as 

dyj 



J(sj,h;hj) 
S(k a ,kj } /j,;jj, s ) 



2ti 
dy 
2tt 
dyk a dy k j 
Att 2 



^B\{x,y^), 



i iyk " ka+iyk J kj S(y ka ,y kj ,ti;n s ). 



At tree-level the position-space jet, beam, and soft functions are given by 

J q{0 \yj,fij) = 1, 



jgi(O) 



X 



y a ,vj 



S qi S(l 



X , 



sShXva,yj,i*s) = i. 



At NLO, the corresponding expressions are 



4:71 

OL s Cf 



7 



2tt 2 



3 

1 + x 2 



+ 3 \n{iyjn 2 e fE ) + 2 In 2 (^j/iV E 
ln(l — x 



2tt 
a s T F 



1 — x 

1 + x 2 



lax + (1 — x) + 



1 — x 

+ \n 2 (iy taf i 2 e~ fE )6(l-x 



(l+x 2 ) 



2tt 



v 1 -x, . 
- ln{iy ta fie lE ){l - 2x + 2x 2 ) 



+ (1 -2x + 2x 2 )(ln 



I — x 



-D + 1 



X 



sSLXyka>ykj,fj-) 



a s C F 
4tx 



4\n 2 (iy ka ^ E ) + 41n 2 (z7/ fcj /2e 7B ) + tt 2 



In arriving at these results in position space we made use of the identities 

-6{z) 



dz e~ lzy 



dz e~ izy 



z 

6(z) In z 



1E\ 



\n(iye 



Iln 2 (^) + ^. 



(B6) 



(B7) 



(B8) 



(B9) 



Appendix C: Renormalization group evolution 



In this section we collect useful formulae that were used in determining the RG evolution 
of the various quantities in the factorization formula given in Eqs. ( 21 ) and ( 22 ) . In particular, 
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we collect formulae for the RG evolution of the hard (H), beam (B), jet (J) and soft (S) 
functions. 

1. Hard function 



The anomalous dimension jh of the hard function is defined by 

^H(Q 2 ,fx) = 7 H#(Q 2 ,/i), 



(CI) 



and can be written as 



I * 



(C2) 



where 7 C is the anomalous dimension of the Wilson coefficient C(Q 2 ,fi) which satisfies 
H(Q 2 ,{i) = \C(Q 2 , fi)\ 2 . The general form of the anomalous dimension 7 C is 

7c = V, ^7cu Sp («,) In -4- + (C3) 



where Sy = 2a^pi ■ pj + iO and <7y = +1 if the momenta Pi and Pj are both incoming or 
outgoing and cr^ = — 1 otherwise. 7 CU s P is related to the cusp anomalous dimension in the 
the fundamental and adjoint representations (a s ) and (a s ) respectively as 



T F (a 

1 cuspl"* 



r A (a 



7cusp \(%s ) ■ 



(C4) 

U F <^a 

The cusp and non-cusp anomalous dimensions and the beta function have expansions in a s 
given by 



7cusp [O^s 



n=0 



4tt 



n=0 



oo 

■2a, £ ( 



n=0 



47T 



(C5) 



For NNLL resummation we need 7 CU s P [65, 66J, 7* [67J, and /? [SHI EH] to 3-loops, 2-loops, and 
3-loops respectively along with NLO PDFs. The 1-loop, 2-loop, and 3-loop cusp anomalous 
dimension coefficients are given by 



^cusp 


= 4, 




^cusp 


= 4 


/67 

At 


^cusp 


= 4 





)C A -—T F n f 

/2i:> 134vr 2 llvr 4 22 
~27~~ + ~45~ + T^ 3 



+ C F T F rif 



55 

~3~ 



I6C3 ) -=^T 2 F n 2 f 



16, 
27' 



M 27 27 3 



(C6) 
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and the beta function coefficients up to 3-loops are given by 

Po 
Pi 

$2 



Y ~~ 3 ' fl 

^C\-™C A T f n f -AC F T F n h 

+ T f n fi 2C F - *fc F C A - ^ CD + Ty/^ + fc A ). 



(C7) 



We define two useful quantities S(/if,Hi) and A(f/,f, m) are needed for the evolution of the 
hard, jet, beam, and soft functions as 



A(ji f ,fii) 



da ^ 



da' 



as{ vf) da 



as{m) P\- a \ 



7cusp ]P^\ 



(C8) 



The expansion of these quantities in a s up to terms needed for NNLL resummation are given 
by 



cusp 

7o 



^ \ lnr ) + 



Afil I a s {ni) \ r 



7l cusp p 

cusp 

7o 



— (1 — r + lnr) + —!r ln 2 r 
Po ' <Zpo 



+ 



aJjM 
Air 



Pili Pi 



Polo nsp Po 



[1 — r + r lnr) + 



(I ~> 



r) lnr 



cusp cusp 



PI _ P2 _ PllT- , 12 

pi Po Poir 



cusp 1 cusp 



7o 



(1 -rf 



and 



A{nf,Hi 



cusp 

7o 
2/?o 



logr + 



otsita) ( li 



4tt V7o 



cusp 



cusp 



Pi 

Po 



^) P2 Pi fi c r 



16tt 2 



7o cusp A) Po \% sp 



(r-1) 

Pi 
Po 



r 2 -l 



The solution to the RG equation in Eq.(Cl) gives the evolution factor 



(C9) 



(CIO) 



H(Q 2 ,n,fi H ) = U H (Q 2 ,n,n H )H(Q 2 ,n H ), 
U h {Q 2 ,H,Lih) = exp 4:C F S(fjL,fiH) -2A H (n,ii H ) 



yo 2 ) 



(Cll) 
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where 



A H (fif,fii 



a '^' da 



IhVA- 



(C12) 



and the 1-loop and 2-loop non-cusp anomalous dimensions for quark fields are given by 

Qfil 1 1 7T 2 9fiD 47T 2 

7^ = C7|,(-3 + 4tt 2 - 48C 3 ) + C F C A {~ ~^- + 52C 3 ) + C^n^ + 



(C13) 



The expansion of Anif^f, fii) is given by replacing 7oJ p — >■ Jh Hi * n Eq.(ClO). 



2. Beam, jet, and soft functions 

The RG equations for the beam, jet, and soft functions are given by the convolution 
equations 

fM-^B q A (x,t,fM) = J dtf J B (t-t',fM) B%(x,t',ii), 

fi-^J(s,fi) = J ds' 7j(s- s',//)J(s',//), 
^■^S{k a ,kj,ii) = J " dk! a J dk'j-f S (k a - k' a ,kj - k'j,fx)S(k' a ,k'j,fx), (C14) 
where the anomalous dimension for the soft function 75 takes the separable form 

7s(ka,kj,fj) = S(k a )-f S (kj,iJ l ) + 5(kj)js(k a ,iJ l ). (C15) 

The anomalous dimensions for the jet, beam, and soft functions have the general form 

1 (n 2 6(s) 



7j(s,/i) = -2C F 7 cusp (a s ) 



7s(t,/u) = -2C F 7 cusp (a s ) — 



+ j q B (a s )5(t), 



7 5 (/c,/i) = 2C F 7 cusp (a s ) +J s (a s )5(k), 



(C16) 



where we have defined the scale jl = [i\/s a j in the soft function anomalous dimension. 
It is often simpler to work in the Fourier transformed space of the beam, jet, and soft 



functions. For example, the factorization formula in Eq.(41) is expressed in terms of the 
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Fourier transformed quantities. The beam, jet, and soft functions and their position space 
analogs are related by 

dy_ 
dyj 



B q A (x,t,n) 
J(sj,h;hj) 
S(k a ,kj,n;fi s ) 



j ty B\(x,y^\ 



2tt 
dyk a dyi 



eWJ^Kfij), 



4tt 2 



±jy ka ^y k jKj S (y ka ,y kj ,n;ns), 



Going into position space, the RG equations take the simpler form 



V-^B g A (x,y,fj,) = -f B (y,fi)B A (x,y,fi), 

Hj^S{y a ,yj) = is(y a ,v) + is(yj,ii) S(y a ,yj), 

where the position space anomalous dimension is defined as 

lB(y,fJ>) = / dte- ity lB {t,n), 



lj(y,v) = J dse iys -fj(s,fi), 
ls{v,») = I dke- iky ls {k,n). 



(C17) 



(C18) 



(C19) 

These position space anomalous dimensions take the general form 

1b(v,v) = 2C F 7 cusp (a s )ln(i?//i 2 e 7B ) +7 9 (a s ), 
1j(y,n) = 2C F 'j cusp (a s )\n(iyn 2 e' yE ) +7 <? (a s ), 

ls (y,n) = -2C Flcusp (a s )\n(iyfie^)+ r s (a s ). (C20) 
The evolution equations in position space are given in terms of the evolution factors C/j as 

B q A (x,y,n;ii B ) = U B (y, V, Hb)B a (x, y, fi B ), 
J(y,H;Hj) = Uj(y,n,n B )J(y,Hj), 

S(y a ,yj,v;vs) = U s (y a ,yj, fj>,fj>s)S(y a ,yj, vs), (C21) 

are are given by 



U B (y ta ,ii f ,Hi) = exp -AC F S{fif,fii)-A B (nf,Hi) (iytalfie"*) 



Uj(y,nj,Hi) = exp - AC F S(n f , m) - Aj(fj, f , m 



%y[i\e m 



Us(y a ,yj,n,iis) = y a yj{i^s€ nE \/§7j) 



2C F A{^ S ) 



exp 



4C F S(n,ns)-A s (n,ns) 



(C22) 
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where we have defined the quantities 



A s {n,ns) 



and 



da 



7b H» 



a.( W ) £H 
•/a.( W ) PFJ 

/ 7?uT 7s ^' 



(C23) 



Tb = Tj, 



9 <7 Q 

Is = -Tj -Tb~ Th- 



(C24) 



The a s expansion of 7} is given by 



oo 



n+l 



7j n ' 



n=0 



and the terms needed for NNLL resummation are 



(C25) 
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.121 2tt 2 . 



— - mCz)C A + (3 - 4tt 2 + 48( 3 )C F + (— + — )/3 
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